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Abstract 

This article deals with the numerical calculation of eigenvalues of per- 
turbed periodic Schrodinger operators located in spectral gaps. Such op- 
erators are encountered in the modeling of the electronic structure of crys- 
tals with local defects, and of photonic crystals. The usual finite element 
Galerkin approximation is known to give rise to spectral pollution. In this 
article, we give a precise description of the corresponding spurious states. 
We then prove that the supercell model does not produce spectral pollu- 
tion. Lastly, we extend results by Lewin and Sere on some no-pollution 
criteria. In particular, we prove that using approximate spectral projec- 
tors enables one to eliminate spectral pollution in a given spectral gap of 
the reference periodic Schrodinger operator. 

1 Introduction 

Periodic Schrodinger operators are encountered in the modeling of the electronic 
structure of crystals, as well as the study of photonic crystals. They are self- 
adjoint operators on //^(M'^) with domain H^(R'^) of the form 

^pcr = - A + Vpor, 

where A is the Laplace operator and V^er a 7^-periodic function of Lj'^j,(R'') {TZ 
being a periodic lattice of E"*), with p = 2 if d < 3, p > 2 for c? = 4 and p = d/2 
for d > 5. 

Such operators describe perfect crystals, by contrast with real crystals, in 
which the underlying periodic structure is perturbed by the presence of local 
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or extended defects. In solid state physics, local defects are due to impurities, 
vacancies, or interstitial atoms, while extended defects correspond to disloca- 
tions or grain boundaries. The properties of the crystal can be dramatically 
affected by the presence of defects. In this article, we consider the case of a 
d-dimensional crystal with a single local defect, whose properties are encoded 
in the perturbed periodic Schrodinger operator 

i? = + = -A + V^pcr + W', WeL'^im.'^), W{x) 0. (1) 

|a:|— foo 

Note that we do not assume here that W is compactly supported. This allows 
us in particular to handle the mean- field model considered in |6]. In the latter 
model, d — 3 and the self-consistent potential W generated by the defect is of 
the form W — p * \ ■ with p e L^(R^) n C, C denoting the Coulomb space. 
Such potentials are continuous and vanish at infinity, but are not compactly 
supported in general. 

Computing the spectrum of the operator H is a key step to understand the 
properties of the system. It is well known that the self-adjoint operator H^^j. 
is bounded from below on L^(R'^), and that the spectrum a{Hp^j.) of i?por is 
purely absolutely continuous, and composed of a finite or countable number of 
closed intervals of K [16] . The open interval laying between two such closed in- 
tervals is called a spectral gap. The multiplication operator W being a compact 
perturbation of H^^j., it follows from Weyl's theorem jl6j that H is self-adjoint 
on L^(R'^) with domain iJ^(R''), and that H and H^^j. have the same essential 
spectrum: 

Contrarily to H^^^., which has no discrete spectrum, H may possess discrete 
eigenvalues. While the discrete eigenvalues located below the minimum of 
(Tcss(-ff) are easily obtained by standard variational approximations (in virtue of 
the Rayleigh-Ritz theorem tl6j), it is more difficult to compute numerically the 
discrete eigenvalues located in spectral gaps, for spectral pollution may occur [5]. 

In Section [21 we recall that the usual finite element Galerkin approximation 
may give rise to spectral pollution [5 , and give a precise description of the 
corresponding spurious states. In Section [3l we show that the supercell model 
does not produce spectral pollution. Lastly, we extend in Section 2] results by 
Lewin and Sere [14| on some no-pollution criteria, which guarantee in particular 
that the numerical method introduced in [6], involving approximate spectral 
projectors, and is spectral pollution free. 
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2 Galerkin approximation 



The discrete eigenvalues of H and the associated eigenvectors can be obtained 
by solving the variational problem 

find {i^, A) e W^iW^) X R such that 

where (•, ■)]^2 is the scalar product of L'^iW'-) and a the bilinear form associated 
with H: 

A sequence (X„)„gN of finite dimensional subspaces of H^{R.'^) being given, we 
consider for all n E N, the self-adjoint operator H\x^ ■ Xn — > X„ defined by 

V(?A„,0„) e Xn X Xn, {H\x„'(pn,4'n)L^ = a{lpn, (l>n) ■ 

The so-called Galerkin method consists in approximating the spectrum of the 
operator H by the eigenvalues of the discretized operators -ff |x„ for n large 
enough, the latter being obtained by solving the variational problem 

find (-0„, A„) e Xn X M such that , , 

V(/)„ e Xn, a{^n,4'n) = A„ (?A„ , (/)„ )i2 . 

According to the Rayleigh-Ritz theorem [TFJ , under the natural assumption that 
the sequence {Xn)nen satisfies 

'icp^H^W'), inf U~cpn\\m^Q, (3) 

the Galerkin method allows to compute the eigenmodes of H associated with 
the discrete eigenvalues located below the bottom of the essential spectrum. It 
is also known (see e.g. [S] for details) that, as H is bounded below, ([3]) implies 

cr(ff) C liminfCT(ff|x„), (4) 

n— >oo 

where the right-hand side is the limit inferior of the sets a {H\x„), that is the 
set of the complex numbers A such that there exists a sequence (A„)„gN, with 
Xn G '^{H\xn) for each n G N, converging toward A. In particular, any discrete 
eigenvalue A of the operator H is well-approximated by a sequence of eigenvalues 
of the discretized operators H\x„- On the other hand, ([3]) is not strong enough 
an assumption to prevent spectral pollution. Some sequences of eigenvalues of 
a{H\x^) may indeed converge to a real number which does not belong to the 
spectrum of H: 

limsupcr (iJ|x„) ^ cr(^^^) in general, (5) 

n— ^oo 

where the limit superior of the sets a (iJ|x„) is the set of the complex numbers A 
such that there exists a subsequence {a{H\x„^ ))fcGN of (cr(iJ|x„))rieN for which 

Vfc e N, 3\nk e (t{H\x^ ) and hm A„^^ A. 
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Spectral pollution has been observed in many situations in physics and me- 
chanics, and this phenomenon is now well-documented (see e.g. [9] and ref- 
erences therein). In [5 , Boulton and Levitin report numerical simulations on 
perturbed periodic Schrodinger operators showing that "the natural approach of 
truncating M"* to a large compact domain and applying the projection method to 
the corresponding Dirichlet problem is prone to spectral pollution!^ . Truncating 

indeed seems reasonable since it is known that the bound states of H decay 
exponentially fast at infinity [15]. The following result provides details on the 
behavior of the spurious modes when the approximation space is constructed 
using the finite element method. 



Proposition 2.1. Let (7^i°°)rieN be a sequence of uniformly regular meshes of 
M.'^, invariant with respect to the translations of the lattice TZ, and such that 
hn ■— max/f g7-^ diam(i4r) ^-n^oo 0. Let (n„)„gN be an increasing sequence of 
closed convex sets ofW^ converging to M'^, Tn ■= {K G \ K C Vln} and Xn 
the finite- dimensional approximation space of H}^(p.n) ^ Ll^iW^) obtained with 
Tn and Pm finite elements (m £ N*^. Let A £ limsup„^oo tr (i/|x„) \ (^{H) 
and (ipUk^Kk) G Xn^ xR be such that H\x„^'4'nk = K^'tpuk, W'^PuJIl^ = 1 and 
limfe_^oo '^Kfc — A. Then, the sequence {ipnkikeN, considered as a sequence of 
functions ofH^{W'^), converges to weakly in H^{R'^) and strongly in L'^^^{W^), 
with q = oo if d = 1, q < oo if d — 2 and q < 2d/ {d — 2) if d > 3, in the sense 
that 

Vi^cR'', K compact, / liJ^J"^ — > 0, 

Jk k^co 

and it holds 

Ve > 0, 3R>0 s. t. liminf / |V^„j2 > 1 - e. (6) 



The latter result shows that the mass of the spurious states concentrates on 
the boundary of the simulation domain 57„j. . 

This phenomenon is clearly observed on the two dimensional numerical sim- 
ulations reported below, which have been performed with the finite element 
software FreeFem-|— I- [11], with Vpcv{x,y) = cos(x) -I- 3sin(2(x -I- y) -I- 1) and 
W{x,y) — — (a;-|-2)^(2y — l)^exp(— (x^-l-y^)). We have checked numerically, us- 
ing the Bloch decomposition method, that there is a gap (a, (3), with a ~ —0.341 



and (3 ~ 0.016, between the first and second bands of iJpor 



-A -I- Vpcr. We 



have also checked numerically, using the pollution free supercell method (see 
Theorem 13. 1 1 below) . that H = Hp^^. + W has exactly one eigenvalue in the gap 
(a,/3) approximatively equal to —0.105. Our simulations have been performed 
with a sequence of Pi-finite element approximation spaces (Xn)4o<n<iooi where 
for each 40 < n < 100, 



rn„ rUn 

-An , 47r 

n n ■ 



with 



40 



20 
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• 7^j°° is a uniform 27rZ^-periodic mesh of consisting of isometrical 
isoceles rectangular triangles per unit cell. 

The spectra oi H\x„ in the gap (a, /3) for 40 < n < 100 are displayed on Fig. 1. 
We clearly see that all these operators have an eigenvalue close to —0.1, which 
is an approximation of a true eigenvalue of H. The corresponding eigenfunction 
for n = 88 (blue circle on Fig. 1) is displayed on Fig. [5] (top); as expected, 
it is localized in the vicinity of the defect. On the other hand, most of these 
discretized operators have several eigenvalues in the range (a,/3), which cannot 
be associated with an eigenvalue of H, and can be interpreted as spurious modes. 
The eigenfunction of H\x„ close to —0.290, obtained for n = 88 (blue square on 
Fig. 1), is displayed on Fig. [5] (bottom); in agreement with the analysis carried 
out in Proposition 12.11 it is localized in the vicinity of the boundary of the 
computational domain. 
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Figure 1: Spectrum of H\x„ in the gap {a, (3) for 40 < n < 100 



Remark 2.1. Using the results in \19i/ . it is possible to characterize the spurious 
states generated by finite element discretizations of one- dimensional perturbed 
Schrddinger operators: for TZ ~ bZ, and r2„ = [— (n + t)b, (n + t)b], the spuri- 
ous eigenvalues are the discrete eigenvalues in [min(cr(_ffpj,j.)), +oo) \ cr(-ff) of 
the operators H+{t) and {t) on L'^{^+) with domains H'^{^+) ni7(J(K+), 

respectively defined by H^{t) — ~-r^ + V^cr(a; ± tb). Besides, the spurious 
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Figure 2: A true eigenfunction, localized close to the defect (top), and a "spu- 
rious" eigenfunction, localized close to the boundary (bottom). 

eigenvectors of H\x„ converge (in some sense, and up to translation) to the 
discrete eigenvectors of H^{t). As 

I y a(i?±(t)) j n[min(a(i/0„)),+^) = [min(a(i/0„)),+oo), 

any A G [niin(cr(iJp(,j.)), +cxd) \ a{H) is a spurious eigenvalue, in the sense that 
there exists an increasing sequence (rin)neN of closed intervals o/R converging 
to M such that 

A G lim inf (7(i/|x„ )■ 

We refer to UOf for a proof and a numerical illustration of this result. The proof 
of similar results for d > 2 is work in progress. 
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Proof of Provosition \2.1\ Wc first notice that, since H — -^A+j (—A + 2Vpor) + 
W , with W bounded in L°°(K'^) and —A + 2V^cr bounded below, there exists a 
constant C G K+ such that 

V^Gi?i(M'^), a{i^,i>)>\\\V^\\l,-C\ml2. (7) 

As 

Vfc G N, IIV'riJlL^ = 1 and a(V'«fc,'0nJ = A„j^ — > A, 

A:— >oo 

we infer from JT]) that the sequence {'iljnk)km is bounded in ff^(M''). It therefore 
converges, up to extraction, to some function G /f^(M''), weakly in H^{W^)^ 
and strongly in L'^^J^'^) with g = c)oifd=l,(7<ooifc?=2 and q < 2d/(d—2) if 
d > 3. It is easy to deduce from ^ and the continuity of a on x i/i(R'') 

that satisfies H(j) = Xcj) and therefore that = since A ^ o'{H) by assumption. 
Consequently, the whole sequence {ipnk)k&i converges to zero weakly in H^{R'^) 
and strongly in i^^^(R'^). 

Let us now prove ([5]) by contradiction. Assume that there exists e > such 
that 

Vi? > 0, liminf / l^PnJ^ < 1 - e. 

'^^oo Jan„^+B{o,R) 

As IlV'nfc ||l2 = 1 for all k, the above inequality also reads 
Vi? > 0, limsup / \tpnj'^ > e, 

fe-i-OO JOH 

where f2^^ = {a; G ri„^ | (i(a;,9il„^) > i?}. We could then extract from {^pnk )keN 
a subsequence, still denoted by (■0„^)fegN, such that there exists an increasing 
sequence {Rnk)kefi of real numbers going to infinity such that 



Vfc G N, 




Let us denote by 

c°{rn = {ve c°(R'^) I yK G TT, v\k g p™} 

the set of continuous functions built from and Pm-finite elements, and by 

= c°(r„°°)niJi(R'*). 

The space X!^ is an (infinite dimensional) closed subspace of iJ ^ (R'^) . Obviously 
Xn ^ X!i^. We then introduce a sequence {Xnk)keN of functions of C^(R'') 
such that for all k G N, 

Supp(x„J C Xk = lon andV|a| < (m+1), 1|9"x»J|l~ < 
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for a constant C G independent of k. Let ■i/'nt = PnkiXnk''l'nk)^ where P„^ 
is the interpolation projector on Xn^. For all k gN, \\tpnk\\L^ > ^^^^ and for all 

(a- A„J(^„^,^~) = (a - A„J(x„>„,,'/'^J 

- (a - A„ J (X„, V'n, - J'n, (Xnfc V'n J , 
-{a- A„J(x„,Vn. - PndXn,1pn,),<f>Z) 

- [ (Ax„,Vn.C + 2C Vx„. • V^„J 

= (a - An, ) (iZ-n, , Xn, G^-^ - (Xn. (j^Z ) ) 

-(a - A„J(X„, '(/'„, - PnAXnJ'nJ, cf>ZJ 

- [ ( Ax„, Vn. C + 2C VXn. • V J , 

where we have used that (a - A„J(^„^ , P„^ (Xn^^^J) = since Pn^Xn^^^J e 
. Denoting by 

we end up with 

- A„, ) (V^„, , CJ = (« - ) (Vn. , Xn, 4>Z - iXn, C )) 

- (a - An, ) (Xn, V-n. - -Pn. (Xnfc V'n J , ) 

-/ M/V^n.C- (8) 

Besides, for hn^^ < 1, 

VC e XZ, WXnJZ - ^n.(Xn.OI|ffi < CK.RZUZWh^^ (9) 

for some constant C independent of A: and . To prove the above inequality, 
we notice that for all K G (Xn.CJk e C°°{K), and S'^CJ^ = if 



8 



= m + 1, so that 



^r~~Z \a\—m-\-l 



~Z a — m+l ^ — ' ^ ' 

where we have used inverse inequahties and the assumption that the sequence 
of meshes (7^°°)„gN is uniformly regular, to obtain the last but one inequality. 

Using the boundedness of {i/jnk)keN in H^(M.'^), the properties of Xnk ^nd W, 
and the fact that {tpnk)keN strongly converges to in Lj^q^(M''), we deduce from 
® and (HI) that 



where the sequence of positive real numbers {rink)keN goes to zero when k goes 
to infinity. 

We can now use Bloch theory (see e.g. [16]) and expand the functions of as 



where T* is the first Brillouin zone of the perfect crystal, and where for all 

For each g G F*, the function {4>Zk)q belongs to the complex Hilbert space 
Ll{T) := {w(x)e*«-=', v e « 7e-periodic} , 



where F denotes the Wigner-Seitz cell of the lattice TZ (notice that the functions 
{(pZ )q are complex- valued) . Recall that if 7?, = hU^ (cubic lattice of parameter 
6 > 0), then F = {~h/2,h/2Y and F* = {-n jh^ix jh^ . The mesh TZ"^ being 
invariant with respect to the translations of the lattice ?<!., it holds in fact 

(0,TJ.eC^°('7;.T)ni^'(r). 
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We thus have for all (j)^^ £ , 

(a° - A„J(V^„,,C) = i «" A„)((4j„ {cl^Z)i) dq. 

where 
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p. 



= VV'J • + VrV',*'/',- (10) 

Let {en,i,q,en,i,q)i<i<N„, £n,i,q < £71,2,9 < ' ' ' < ^n,N^,q, be an (r)-orthonormal 
basis of eigenmodes of a" in C°(7^j°°) fl ig(r). Expanding ('0„fc)g in the basis 
(enfc,/,g)i<i<Ar„ , we get 



i=i 



Choosing 0^ such that 



J = l 



we obtain ||(/'^J|l2 = H^-nJU^ and 



It is easy to check that liminf max |e„j,.j_g — A„^| — Q :— dist(A, (T(-ff" > 0. 
Hence, _ 

liminf(aO-A„J(V'„,,C) >Ce- 

Besides, 

IICIU^ = II^^"JU^ and a°(0-,CJ=«°(V^n.,V^nJ, 
which implies that the sequence ((?!)5^)fegN is bounded in H^{W'-). Consequently, 

< Ce < liminf(a" - A„J(V^„,,0-) < liminf ,7„J|<^- = 0. 

fe— S-oo fc— Kx; 

We reach a contradiction. □ 

A careful look on the above proof shows that the assumptions in Propo- 
sition 12.11 can be weakened: in particular, the mesh Tn can be refined in the 
regions where \ W\ is large, and coarsened in the vicinity of the boundary of r2„ 
(see [To] for a more precise statement). 
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3 Supercell method 



In solid state physics and materials science, the current state-of-the-art tech- 
nique to compute the discrete eigenvalues of a pcrtiirbed periodic Schrodinger 
operator in spectral gaps is the supercell method. Let TZ be the periodic lattice 
of the host crystal and F its Wigncr-Scitz cell. In the case of a cubic lattice of 
paramater > 0, we have TZ = bU^ and F = (—6/2, 6/2]''. The supercell method 
consists in solving the spectral problem 

J find (V'L.Af, Al, at) G Xi^AT X M such that 

\ ^<t>L,N & Xl,N, aL{i'L,N,4'L,N) = '>^L,n{'>I^L,nAl,n)lI^^{Tl)i 

where F^ = LF (with L € N*) is the supercell, 

'^^per(rL) = {uL € L^^iM.'^) \ ul L7^-periodic} , 

and X^^jv is a finite dimensional subspace of 

We denote by Hl,n = HL\xL^fi, where Hl is the unique self-adjoint operator 
on Lpg^(Fi) associated with the quadratic form oi,. It then holds that D{Hl) = 

and 

V(/)L,Ar G Xl,N, Hl,n4>L,N = -^<pL,N + ^Xl.n ((^pcr + Wl)4>L,n) , 

where Wl G L'^^^iTj^) denotes the L7?,-pcriodic extension of iy|ri, and 11^^, „ is 
the orthogonal projector of L'^^^{Tl) on Xl,n for the L-^^^iVi,) inner product. 

Again for the sake of clarity, we restrict ourselves to cubic lattices iJZ = b'E'^) 
and to the most popular discretization method for supercell model, namely the 
Fourier (also called planewave) method. We therefore consider approximation 
spaces of the form 



^L,N = { ^ CkeL,k I Vfc, C-k = C*k 

^fe£27r(6L)-iZ<i I |fe|<27r(bZ,)-iJV 

where clA^) = \T L\^'^/^e'^-=' . 

From the classical Jackson inequality for Fourier truncation, we deduce by 
scaling the following property of the discretization spaces Xl^n: for all real 
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numbers r and s such that < r < s, there exists a constant C > such that 
for all L G N* and all € H^c^{Tl), 

Ul - Ilx,,MH;^Ar.) < ^ \\M\H;^Ar.y (12) 

Our analysis of the supercell method requires some assumption on the potential 
Vpci- We define the functional space A^pcr(r) as 

{||t;w||L2 (pj^) ] 
v e L (T) I ||t;||Aip„(r) := sup sup n <oo \ . 

LeN-„,eHi^,,(rx,)\{o} iMlH^^A^r.) J 

It is quite standard to prove that A^por(r) is a normed space and that the space 
of the 72.-periodic functions of class C°° is dense in A^por(r). We denote the 
72.-periodic Lorentz spaces [4] by LP'^j.{T). 

Proposition 3.1. The following embeddings are continuous: 
ford-1, Ll^^iT)^ Mpor{r), 

ford = 2, L%^ir)^Mpcrir), 

ford = 3, L%^{r)^Mp,A^). 

Proof. We only prove the result for d = 3; the other two embeddings are ob- 
tained by similar arguments. Let us first recall that the Lorentz space L^'°°{T) 
is a L^-multiplier of L^'^(r) (this can be seen by combining results on convolu- 
tion multiplier spaces [2] and continuity properties of the Fourier transform on 
Lorentz spaces [4l), in the sense that 

3Ci £ R+ I V/ G L^''^{T), Vg G L^^'^iT), \\fg\\mr) < Ci||/|U3,=o(r)||g||L«.^(r). 
Besides, the embedding of H^{r) into i^'^(r) is continuous (see for instance) 
3C2 G M+ I Vg G H\T), Wgh^.^^r) < C2\\g\\mir)- (13) 
Let V G ipef (r). Denoting by II := 7^ n {-Lb/2, Lb/2f, we have, for ah 
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li?er(ri,) 

^ j \v{x)w{x + R)\Ux= J2 \\vw{. + Rm 

< C\ Mli3,oo(r)||«^(. + i?)||ia,.(r) 

< C^flklli3.oo(r) E ll^(- + ^)lli^.^(r) 

|3,oo(r) 

-RGli, 

< C^?t^2lklli3.~(r) E /(K^ + ^^)l' + |V«;(x + i?)|2) da; 

< C'i2C'|||z;||i3,o.(r) / (Mx)p + |Vu;(x)p) da; 

< C'2c|||^;||i3.o„(r)|k||ii^,(r,)- 

Therefore, w G Alper(r) and ||t^||Mp„(r) < C'iC2||w||L3,oo(r). □ 

Remark 3.1. In dimension 3, the TZ-periodic Coulomb kernel G\ defined by 

-AGi = 47r ( E -^fl - |rr' ) , mm dix) = 0, 

is in L^'^{T), hence in AdpeA^)- The functional setting we have introduced 
therefore allows us to deal with the electronic structure of crystals containing 
point-like nuclei. 

Theorem 3.1. Assume that Vpor G A^por(r). Then 

Urn a{HL,N) = a{H). 

Proof. Let us first establish that 

o-(H) C liminf a(HL n)- 

Af,L->oo I Af/L^oo 

Nl 

Let A G a-{H) and {Nl)l&N' be a sequence of integers such that — — — > oo. 
Let e > and V e C^°°(M'*) be such that ||V'||z,2 = 1 and \\{H - A)V'||z,2 < e. 
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We denote by '0L the L7?.-periodic extension of i/'Iti,- Since ip is compactly 
supported, there exists Lq G N* such that for aU L > Lo, S'wpp{ip) C Tl. 
Consequently, for aU L > Lq, 

IIV'illLj^irz.) = 1 and \\{Hl - A)V'l ||L2_^^(r,,) < £• 
Let iPl,Nl Ilxf, Ni^i^L- We are going to prove that 

\\{Hl - X)^L - {Hl.n, - A) V^L.AT, 11^. (p . 0. (14) 

First, we infer from ([T^ and the density of -ffo(f^) in L^(f2) for any bounded 
domain Q, of R*^ , that 



e L 



2lmd\ 



(1 - nx_t.„ (Fx,) ^ — > 



0, 



where L^iW^) denotes the space of the square integrable functions on with 
compact supports, and where 4>l is the LT^-periodic extension of 4'\tl- tp, 
Alp, Vpertp and Wip are square integrable, with compact supports, we therefore 
have for all L > Lq, 



L 



0, 



\-A^L + A^L.N, l|L^„(r,) = II (l - nx,,„, ) (-A^) 
IWLi^L - Tlx,,., iWMhi^^ir,) = II (l - nx,. 



0, 



0, 



0. 



We infer from the last two convergence results that, on the one hand, 



< 



< 



WL4'L-Tlx,,.,iWL^L] 

WLi>L-Tlx,,,,{WL^L] 
0, 



+ 



n 



Xl, 



{WL{i^L-i^L,N.)) 



+ 



L—^oo 

and that, on the other hand, 

IIVpcrV'L - TlxL,N,{ypcr1pL,NL)\\Ll^A'ri^) 



VperIpL - Tlx L.N, (Vpcr-^L) 
Vpci-IpL 



TlXr.,Nr (^PorlV'i " V'L.JVi,)) 



L^OO 



1 - TIXl,m, 
0. 



+ ||Vpor|Up„(r)||'0L - V'L.ATi.llffi (Fx,) 
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Collecting the above results, we obtain Thus, for L large enough, 

||(ifL,jV, -A)VL,ivJ|L^„(r,) <2£. 

As IIV'L.ATi, ||l2 (Pj^) = 1 for all L > Lq, we infer that for L large enough, 
dist{X, a{HL,NL)) < 2e, so that A G hminf cr(iJL.Ar^). 

Let us now prove that 

limsup (j{Hn,l) C (t{H). 

N,L^oc \ N/L-^oo 

We argue by contradiction, assuming that there exists A £ K \ (j{H) and a 
sequence (Lk, Nk)k(ETi with Lk — J> oo, Nk oo, Nk/Lk oo, such that 

A:— >oo k-^oo A:— J-oo 

for each fc, there exists {'ipL^.N^j ^l^.n^) 6 Xl^^^^ x R satisfying 

I IIV-Lk.JvJlL^fr^j = 1, 

and lim \Lk,Nk = A. Each function ipLk,Nk is then solution to the PDE 

- \^i'Lk,Nu +nxz,,,„, ((Vr + W^iJV'Lfc,ArJ = Ai, ,JV,'(/'i.,Wfc • (15) 

Reasoning as in the proof of Proposition [^?T1 it can be checked that the sequence 
(IIV'Lfc,wJ|ffi„(ri.j)fceN is bounded, and that 

V'L^Ar.— >0 inLL(K')- (16) 

For all A;, we consider a cut-off function Xk ^ C^(R^) such that < Xfc ^ 1 
on Xk = 1 on r^,, Supp(xfc) C (L^ + i[/')r, ||Vxfc||L- < CLl^'\ and 
||Axfc||ioo < CLf,^ for some constant C e M+ independent of fc. We then set 
i^k = Xk^L^.N,. It holds ^k e i?2(R'^), 1 < W^kh^ < 2''/^ and 

-^AxfeV'L,,Wfc + (AL^.ATfe - A)^fe. (17) 
As {XLk,Nk)kGN converges to A in R and ll-iAfclU^ < 2''/^, we have 
(Al,,jv. - X)i^k strongly in ^^(R'^). 

k-^oo 

— 1/2 

Using the facts that Supp(xfc) C ^Tl^, ||Vxfc[|L~ < CL^ and ||Axfc||L= < 
CL'^'^ for a constant C S R+ independent of fc, and the boundedness of the 
sequence (||^/'Lfc,JvJ|//i„(ri j)feGN, we get 

-Vxfc • VV'L^.A', - \/^Xk^Lk.N, strongly in L''{W'). 
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It also follows from (ITBl) that the sequence \\WLkipLk,Nk\\L^^^{rLj S^es to zero, 
leading to 



Lk.Nk) 



strongly in L^(M''). 



Lastly, 



Xfe(^pcr^L,,iv, -nx,^.„J^pcrV'L..Jvj)^-:;^0 strongly in (M'*) . (18) 

To show the above convergence result, we consider e > and, using the density 
of e.g. lypiir(r) {M^pcr e LS=cr(r) I VW^per € L^er(r)} m 7Wper(r), we Can 



choose some Vpcr G W^i.'^iT) such that ||Vpor - Vpcr||A4per(r) 
deduce from (HI]) that, for all k eN, 



< e. We then 



< 



(v"por?/'L,,Arfc) 



< \\Vpci- - V'pcr||A^p„,(r)II^Lfc,7vJ|Hi^^(ri,J + Tpll^por'0ife,riJ|Hi„(rz,J 



< £|IV'L„ivJ|ifi„(r.J + ^IIV-L.,™. llHi„(r.j(ll VrlU- + l|Vyper||L~). 

Since the sequence ( \\ipLk,Nk\\H^ jTl ) ) bounded, this yields 



VpcrlpLk.Nk ~ ^Xl 



ken 

(VporV'L,,ArJ 



0, 



which implies (|T8] 



Collecting the above convergence results, we obtain that the right-hand side 
of (fT7)) goes to zero strongly in L^(M''). Therefore, {il'k/\\i'k\\L^)k£N is a Weyl 
sequence for A, which contradicts the fact that A ^ a{Hp^^). □ 

A similar result was proved in [18 for compactly supported defects in 2D pho- 
tonic crystals, with Vpci- G L°°(M^) and N = oo. In [7], we prove that the error 
made on the eigenvalues and the associated eigenvectors decays exponentially 
with respect to the size of the supercell. We did not consider here the error due 
to numerical integration. The numerical analysis of the latter is ongoing work 
and will be reported in jlOI. 

Note that, if instead of supercells of the form = LT, i £ N*, we had 
used computational domains of the form F^+f = {L + t)T, t G (0,1), we 
would have observed spectral pollution. As in the case studied in the previ- 
ous section, the spurious eigenvectors concentrate on the boundary c^F^+t. In 
the one-dimensional setting {TZ — hi), and for a fixed value of i, the trans- 
lated spurious modes (j)L,N{- — {L + t)b/2) strongly converge in i7[Q^(R), when 
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L goes to infinity, to the normalized eigenmodes of the dislocation operator 

H{t) = + l^<oV^pcr(a: + th/2) + l^>oVr(a: - tb/2) studied in [l^. We 

refer to [10] for further details. 

4 A no-pollution criterion 

Spectral pollution can be avoided by using e.g. the quadratic projection method, 
introduced in an abstract setting in [17 1 , and applied to the case of perturbed 
periodic Schrodinger operators in 5 . An alternative way to prevent spec- 
tral pollution is to impose constraints on the approximation spaces (X„)„gN. 
Consider a gap (a,/3) C M \ a{H^^^) in the spectrum of and denote by 

P — X(_oo,7] (^^pcr) where 7 = ^^y^ and where X(-oo,7] is the characteristic 
function of the interval (—00,7]. 



Theorem 4.1. Let {Pn)neN be a sequence of linear projectors on L^(M'*) such 
that for all n E N, Ran{Pn) C H^{M.'^), and sup^^pj ||P„||£(2;,2) < 00, and 
{Xn)neN 0, sequence of finite dimensional discretization spaces satisfying ^Bj) 
as well as the following two properties: 

(Al) Vn e N, X„ = X+ © X- with X- C Ran(P„) and X+ C Ran(l - P„); 

(A2) sup — ^ 0. 

Then, 

lim a{H\xJr\{a,l3)=a{H)r\{a,P). 



The above result is an extension, for the specific case of perturbed periodic 
Schrodinger operators, to the results in [TU Theorem 2.6] in the sense that (i) 
the exact spectral projector P is replaced by an approximate projector P„, and 
(ii) the discretization space Xn may consist of functions of 7J^(M'') (the form 
domain oiH), while in 14 , the basis functions are assumed to belong to 7J^(M'') 
(the domain of H). 



Proof. From we already know that cr{H) n {a, (3) C liminf„^oo <^{H\x^) n 
{a, (3). Conversely, let A G {limsup^_^^ a{H\x J n {a, I3))\a-{H), and (V'„j,)fceN 
be a sequence of functions of H^(R'^) such that for all fc G N, tpn^ S Xn^, 
\\^^nk\\L^m^') = 1 and {H\x„^ - >^)ip7ik — ^-fe^oo strongly in L^(R'^). Reasoning 
as in the proof of Proposition 12.11 we obtain that the sequence {'ipnk)keN con- 
verges to 0, weakly in iJ^(K'^), and strongly in Ljq^(R'^). Let us then expand i/'nt 

as Ipn^ = +"0,7^ with {l-Pn,)lpn, G X+, and Ip'^ := Pn^n^ G ^n^ 

and notice that 

(a° - A) « ,0-1- (a° - A) ( , ) = (« - A) (V-n, , <) - /" W^V-n. < • 
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Since V+ = (1 - Pn)i^n, G X„ 



((iJ|x„ -A)V„,,(l-P„)V'n.)L^ 



Besides, as W vanishes at infinity, {i>nk)keti converges to in ^^^^(R'^) and 
suPfegN llV'nfe lli^ < 1 + supfegN \\Pnu Wc^L^) < 00, we also have 



" fc->oo 



Therefore, 



Likewise, 



/ 

(a° - A) « , V+ ) + (a° - A) (V" , V+ ) 



0. 



k—>oo 



0. 



Substracting the second equation from thc^ first one, we obtain 
(a" - A)(V+ , V+ ) - (aO - A)(V- , V-^J 0. 
Now, we notice that 

(a"-A)(V^- ) = (a°-A)(P„,7A„,,P„,V«J 

= (aO-A)(P^„,,PVnJ+2(a°-A)(F^„,,(P„, -P)^„ 
+(aO - A)((P„, - P)V„„ (Pn. - P)V'„J, 

and 

(a°-A)(V'+,0 = (a"-A)((l-P„JV'„,,(l-P„J^„J 
= (a°-A)((l-P)^„,,(l-P)7A„J 

+2{a° - A)((l - P)V«., (P - PnJV'nJ 

+ (aO - A)((P - P„JVn„ (P - PnJV-nJ. 

Besides, there exists i]+,rj- > such that for ah e H^{M.'^), 

(aO-A)((l-P)^,(l-P)^) > r?+||(l-P)^||i.(K.), 



-(a° - A)(PV',PV') > r]-\\Py 



lL2(R<i)- 



Thus, 



(a°-A)(V+,V+)-(«°-A)(V-,CJ > min(^+,,?_)||V'„J|i,(„.) 

+2(a°-A)(Vn„(P-P„JV. 



18 



From assumption {A2) and the boundedness of (V'nfc)fceN in we deduce 

that 

which uTiply that ll^n^llLS — > 0. This contradicts the fact that ll'^n^llLS — 1 

k^oo 

for all keN. □ 



The assumptions made in Theorem 14.11 allow in particular to consider ap- 
proximation spaces built from approximate spectral projectors of -ffpor- 
a matter of illustration, let us consider the case when the approximate spec- 
tral projectors are constructed by means of the finite element method. As 
in Section [21 we consider a sequence (7^°°)„gN of uniformly regular meshes of 
M'*, invariant with respect to the translations of the lattice TZ, and such that 
hn ■= max/f g7-~ diam(X) — >n^oci 0, and denote by the infinite dimen- 
sional closed vector subspace of H^{R'^) built from (7^°°)„gN and Pm-finite el- 
ements. Assume that we want to compute the eigenvalues oi H = H^^j. + W 
located inside the gap (a, (3) between the J"^ and {J+iy^ bands of H^^^. Using 
Bloch theory ^16j, we obtain 

P = Xi-o.n]iKor)=£^P,dQ: 

where Pq is the rank- J orthogonal projector on L'i{T) defined by 
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Pq = \^j,q) {^j,ql^ 



where (ej,q, ej,q)jgN* , < £2,9 < ■ • • , is an L^(r)-orthonormal basis of eigen- 
modes of the quadratic form defined by (jlOp . For n large enough, we introduce 

^" ■= 7 {(^n.j,q\ dq, (19) 

where (e„,j\g, e„j,,)i<j<Ar„ , en,i,q < en,2,q < ■ ■ < <^n,N„,q^ istheL^(r)-orthonormal 
basis of eigenmodes of a° in C°(7^°°) fl L'^{T) already introduced in the proof of 
Proposition 12.11 

We have seen in Section [2] that using approximation spaces of the form 

Xn - {^n e X^ I Supp(^„) C Qn} , 

where (ri„)„gN is an increasing sequence of closed convex sets of M'' converging to 
R'^, leads, in general, to spectral pollution. We now consider the approximation 
spaces 

Xn^X+(B X- where X' = P„X„ and X+ = (1 - P„)^«. (20) 
Note that Xn = Xn + PnXn, so that Xn can be seen as an augmentation of X„. 
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Corollary 4.1. The sequence of approximation spaces (X„)„gN defined by h20\) 
satisfies and it holds 

lim a{H\ ^ ) n (a, /3) = a{H) n [a, (21) 



Proof. As Xn = Xn+PnXn witli (X„)„gN Satisfying ([3]), it is clear that {Xn)nm 
satisfies The sequence {Pn)n&i is a sequence of orthogonal projectors of 
L2(M'^) such that Ran(P„) d X^ C H^iW^). Besides, ||P„||£(l2) = 1 since the 
projector P„ is orthogonal. It follows from the minmax principle [16] and usual 
a priori error estimates for linear elliptic eigenvalue problems [3] that 

sup e„ j q — > a and inf e„ , „ — > (3, 

and that there exists C E M+ such that 



l|t^gl|L2(r) = 1 



We conclude using Theorem 14. II □ 



Let us finally present some numerical simulations illustrating Corollary 14.11 in 
a one-dimensional setting, with Vpci{x) = cos(x) + 3sin(2x + 1) and W{x) = 
— {x + 2)'^e~^ . We focus on the spectral gap (a,/3) located between the first 
and second bands of H^^j. — + Vpcr (corresponding to J = 1). Numerical 

simulations done with the pollution- free supercell model show that a ~ —1.15 
and /? ~ —0.65, and that H has exactly two discrete eigenvalues Ai ~ —1.04 
and A2 — —0.66 in the gap [a, f3). 

The simulations below have been performed with a uniform mesh of R cen- 
tered on 0, consisting of segments of length h = 7r/50, and with $7 = [—L,L], 
for different values of L. The sums over TZ have been truncated using very 
large cut-offs; likewise, the integrals on the Brillouin zone have been computed 
numerically on a very fine uniform integration grid, in order to eliminate the 
so-called fc-point discretization errors. The numerical analysis of the approxima- 
tions resulting from the truncation of the sums over TZ and from the numerical 
integration on F*, is work in progress. 

The spectra of the operators H\x„ (standard finite element discretization 
spaces) and (augmented finite element discretization spaces defined by 

([20]) ) are displayed in Figure ID The variational approximation of H in Xn is 
seen to generate spectral pollution, while, in agreement with Corollarv 14.11 no 
spectral pollution is observed with the discretization spaces X„. 
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16 17 18 19 20 21 22 23 




Figure 3: The spectra of the variational approximations of H for various sizes 
of the simulation domain, obtained with standard finite element discretization 
spaces Xn (top) and with augmented finite element discretization spaces Xn 
defined by ^ (bottom). 
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